Reduction of Josephson critical current in short ballistic SNS weak links 
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We present fully self-consistent calculations of the thermodynamic properties of three-dimensional 
clean SNS Josephson junctions, where S is an s-wave short-coherence-length superconductor and A*' 
is a clean normal metal. The junction is modeled on an infinite cubic lattice such that the transverse 
width of the S is the same as that of the A'^, and its thickness is tuned from the short to long limit. 
Both the reduced order parameter near the 5'A'' boundary and the short coherence length depress the 
critical Josephson current Ic, even in short junctions. This is contrasted with recent measurements 
on SNS junctions finding much smaller IcRn products than expected from the standard (non-self 
consistent and quasiclassical) predictions. We also find unusual current-phase relations, a "phase 
anti-dipole" spatial distribution of the self-consistently determined contribution to the macroscopic 
phase, and an "unexpected" minigap in the local density of states within the A'^ region. 

PACS numbers; 74.50. -f-r, 74.80. Fp 



Over the past decade, both experimental and theoret- 
ical interest in the superconductivity of inhoniogeneous 
systems have been rekindled, thereby leading to a reex- 
amination of eyen well-charted areas from the mesoscopic 
point of view.a For example, the Josephson effect in 
a superconductor-normal-metal-superconductor {SNS) 
weak link was known to be the result of the macroscopic 
condensate wave function leaking from the S into the TV 
region. The induction of such superconducting correla- 
tions in the N, the so-called proximity effect, has been 
given a new real-space interpretation through the relative 
phase-coherence of quasipartides, correlated by Andreev 
reflection at the SN interface. H Moreover, the realization 
of the importance of tracking the phase-coherence of sin- 
gle particle wave- functions in proximity-coupled metals 
of mesoscopic size has also unearthed new phenomena, 
such as quantization of the critical current in ballistic 
mesoscopic short SNS junctions at low enough temper- 
ature. H'Q In short clean junctions, as T ^ 0, the critical 
supercurrent If, = eA/h carried by a single conducting 
channel depends only on the superconducting energy gap 
A as the smallest energy scale A < JJxh = hvp/L (in 
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the long junction limit Ic ^ i?Th is setd by the "ballistic" 
Thouless energy i?Th < A, which is a single quasiparticle 
property determined by the Fermi velocity Vp in the TV 
interlayer of length L). Thus, both mesoscopic and "clas- 
sical" clean point-contact SNS junctions, with ballistic 
transport i > L [i \s the mean free path), are predicted 
to exhibit the same IcRn — ttA/c product at T = 0. 
This has been known for quite some time as the Kulik- 
Omelyanchuk (KO) formula,El where i?7v is the Sharvin 
point contact resistance Rn — hj^e^M of the ballistic iV 
region containing M conducting. channels. „ 

Recent experimental activityQ on highly transparentcl 
ballistic short SNS junctions, where both I^ and Rn 
are independent of the junction length, reveals much 
lower values of IqRn than the KO formula (similarly, 
the critical current steps found in an attemptu to ob- 
serve discretized I^ are much smaller than the predicted 



eA/h). However, a proper interpretation of these results 
demands a clear understanding of the relationship be- 
tween relevant energy and length scales. The criterion 
for the short junction limit A < i?Th introduces a "co- 
herence length" of the junction ^o = ^Vp /i^A\ i.e., the 
maximum KO limit can be expected only for L <^ ^o- 
The relation between ksT and i?Th defines the high- 
{ksT > E'xh) versus low- (/csT < E'xh) temperature lim- 
its, which is equivalently expressed in terms of the junc- 
tion thickness as L > ^^r versus L < ^n, respectively, 
with ^AT — hvp /2'KkBT being the normal metal coher- 
ence length. The S,n sets the scale over which two quasi- 
partides in the N , correlated by Andreev reflection, re- 
tain their relative phase coherence (i.e., superconducting 
correlations imparted on the A'^ region at finite tempera- 
ture decay exponentially with ^^v, while at zero temper- 
ature ^AT — > cx) and the condensate wave funetion decays 
inversely in the distance from the interfacala). There- 
fore, the simple exponential decay of Ic ~ exp(— L/^at) 
appears only in the high-temperature limit, while in the 
opposite low-temperature limit ^^r ceases to be a relevant 
length scale and the decay is slower than exponential. 
The paforementioned experiments on clean SNS junc- 
tionaZl are conducted on Nb/InAs/Nb junctions which are 
tuned to lie in the regime where Cs < ^ < Co ^ Civ 
{S,s = hvp/nA is the bulk superconductor coherence 
length). Thus, the large difference between ^g and ^o 
means that there is a substantial Bsmi velocity mismatch 
(typically an order of magnituddj'u) , which must gener- 
ate normal scattering at the SN interface in addition to 
Andreev reflection. This, together with other possible 
sources ofnSGattering at the SN boundary, like-imperfect 
interfacefflO'EHl or charge accumulation layeralJ (typical 
of Nb/InAs contact), cannot be detected by only observ- 
ing the independence of /g and i?jv on interelectrode sep- 
aration (for intermediateQ L). Nevertheless, this is fre- 
quently the criterion used in experimentaj to ensure that 
the transport is ballistic. Therefore, the ideal maximum 



value for Ic could be achieved only for ^g = ^o 
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FIG. 1. Product of the critical current Ic and the normal 
state resistance Rn as a function of the SNS junction thick- 
ness L. Both the S and A'^ are at half-filling in the bulk. 
The value of IcRn is always below the product of the bulk 
critical current in the S leads and the Sharvin point con- 
tact resistance, /^""'-Rsh = 1.45A/e. The right axis measures 
the ratio of the Kulik-Omelyanchuk formula IcRn = vrA/e 
for the clean superconducting point contact {L — > 0) and 
IcRn of our junctions. The inset shows the decay of the or- 
der parameter at the 5'A'^ interface for / = 0, which reaches 
an asymptotic value of about one-half of the bulk A for 
L ~ 2^s [IcRn / ^{z = SN interface) is virtually constant 
for L < 2^s]- 



and with a perfectly transparent interface, where the 
junction thickness satisfies L <S^ ^s- Even in this case it 
is possible that the current in short junctions is smaller 
than expected due to a depressed value of the order pa- 
rameter on the SN boundary when thc_transverse width 
of the S and N regions are the same.E3 Swh junctions 
cannot be treated by simplified approachegau assuming a 
step function for A{z) because the order papimeter varies 
within the S due to the self-consistency.Ejtj 

Here we undertake an idealized study of different in- 
trinsic properties of three-dimensional SNS junctions 
which can be detrimental to Ic, without invoking any 
sample-fabrication dependent additional scattering at 
the SN interface. Two such effects are known: (i) the 
requirements of self-consistency, which become impor- 
tant for specific junction geometries delineated below, 
depresses the order parameter near the SN boundary 
and therefore the current in short junctions; (ii) a fi- 
nite ratio A/fj, (where /j, is the Fermi energy measured 
from the bottom of the band) generates intrinsic normal 
scattering at the SN boundary (without the presence of 
impurities or barriers at the interface). Therefore, even a 
clean junction (with ballistic transport above Tc) might 
not be in the ballistic limittj below the superconducting 
transition temperature T^ unless the filling is tuned to 
the energy of the transmission resonances. Our princi- 
pal result for the evolution of IcRn as a function of L 
is shown in Fig. [l|. The IcRn drops by about an order 
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FIG. 2. Scaling of the current-phase relation l{<j))/lc with 
the thickness of a clean SNS junction. Note that the phase 
change across the junction (j)c at the critical current Ic ~ I{(i>c) 
varies monotonically with the junction thickness, as shown in 
the inset, and is always far below tt, which is the prediction of 
non-self-Gonsistent calculations in both the shorlJj {L <^ ^o) 
and longQ (L ^ ^o) junction limits at T ^ 0. 



of magnitude at L '^ ^jv , thus showing how the charac- 
teristic voltage can be reduced dramatically in moderate 
length junctions, even in the low-temperature limit (to 
which our junctions belong). The reduction of Ic in our 
short junctions is determined by the depression of the or- 
der parameter in the S, as demonstrated by the inset in 
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Fig. |l| where A(z) at the SN interface decreases asymp- 
totically to a limiting value reached at L ^ 2i^5 with 
IcRn / A.{z = SN interface) being nearly a constant for 
L < 2£^s- For the junctions thicker than 2^5 the decay of 
the critical current Ic ^ 1/L scales asO -Exh, while at non- 
zero temperatures and for long enough junctions L > ^^r 
it changes into a simple exponential decay. Thus, in the 
general case ^s < £,n, Ic can be independent of L only for 
2^5 < L < £.0, as observed in the experiments. However, 
such thickness-independent Ic can be substantially below 
MeA/h, as defined by the inevitable {vp ^ Vp) interface 
scattering and/or reduced A, with its lowest value being 
set at L ~ 2^5 by the "inverse proximity effect" on the 
S side of a SN structure. We believe that ballistic be- 
havior could be found in our junctions at even lower T, 
where £_n ^ 2^5, but such calculations are technically 
more involved at present. 

The SNS Josephson junction is modeled by a Hubbard 
Hamiltonian 



H=-J2t, 



cLc, 



10 a 



X]t/J4|Qt--) \c\l<^^i-7; 



(1) 



on a simple cubic lattice (with lattice constant a). Here 
c\^ (cia-) creates (destroys) an electron of spin a at site 
i, tij = t (the energy unit) is the hopping integral be- 
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FIG. 3. Scaling of a spatial distribution of the phase devi- 
ation S(j){z) within the self-consistently modeled part of the 
clean SNS junction. The total phase change across the junc- 
tion (f) is the sum of the bulk phase gradient x L and the change 
in 5(l){z) along the A*' interlayer, Eq. (y): 00 at some small su- 
percurrent and <j)c at the critical junction current Ic- At large 
enough junction thickness L [panel (a)] the shape of 5(l>{z) is 
just rescaled by the increase of the Josephson current, while 
for smaller L the shape changes abruptly upon approaching 
Ic [panels (b) and (c)]. 



tween nearest neighbor sites i and *^ The SNS struc- 
ture is comprised of stacked planealJ where Ui < is 
the attractive interaction for sites within the supercon- 
ducting planes (inside the N region Ui = 0). In the 
Hartree-Fock approximation (HFA), this leads to a BCS 
mean-field superconductivity in the S leads, where for 
t/j = -2 and half-filling we get A = 0.197t (T^ = O.lli) 
and fs — hvp/nA ~ 4a. The lattice Hamiltonian (|l|) 
of the inhomogeneous SNS system is solved by comput- 
ing a Nambu-Gor'kov matrix Green function. The off- 
diagonal block of this matrix is the anomalous average 
which quantifies the establishment of superconducting 
correlations in either the S [A(z) = \U{z)\F{z), where 
F{zi,Zi,T = 0+) is the pair-field amplitude] or the N 
region. For the local interaction treated in the HFA, 
computation of the Green function reduces to invert- 
ing an infinite block tridiagonal Hamiltonian matrix in 
real space. The Green functions are thereby expressed 
through a matrix continued fraction (technical details 
are given elsewherellZl) . The final solution is fully self- 
consistent in the order parameter |A(z)|e*'^^^^ inside the 
part of the junction comprised of the N region and the 
first 30 planes inside the superconducting leads on each 
side of the N interlayer. Our Hamiltonian formulation of 
the problem and its solution by this Green function tech- 
nique is equivalent to solving a discretized version of the 




FIG. 4. Local density of states at the central plane of 
a clean 3D SNS junction (composed of 10 normal planes, 
L — 10a) for different supercurrent flows: 7 = 0, 7 = 7c/4, 
and 7 = 7c/2. The peaks correspond to ABS confined within 
the A'^ region at energies E < A. For 7 7^ 0, the degeneracy 
of right- and left-moving electrons is lifted by a Doppler shift, 
|g rise to the Josephson DC current (or at least part of 
). The minigap around the Fermi energy lj = in the A'^ 
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region appears to be a result of a finite A/fi 
normal scattering at the SN interfaces. 
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Bogoliubov-de GennealS ffldG) equations formulated in 
terms of Green functions Jlj but in a fully self-consistent 
manner — by determining the off-diagonal pairing poten- 
tial A(z) in the BdG Hamiltoniant3 after each, iteration 
until convergence is achieved. The tight-binding descrip- 
tion of the electronic states also allows us to include an 
arbitmry band structure (which is rarely taken into ac- 
countcj), or more complicated pairing symmetries. The 
calculation is performed at T = 0.09Tc where ^n = 40a, 
which is a low-temperature limit for almost all of our 
junction thicknesses. 

This technique is different from the quasiclassical use 
of a coarse-grained microscopic Gor'kov Green function, 
through either the Eilenberger equations (clean limit) or 
Usadel equations (dirty limitI,El or non-self-consistent 
solutions of the BdG equationsH which are applicable only 
for special geometries where the left and right S leads can 
be characterized by a constant phase 0l and (Jjr, respec- 
tively. This neglects the phase gradient {d(j)/dz) hulk in- 
side the S, thereby violating current conservation. Such 
an assumption is justified when the critical current of 
the junction is limited by, e.g., a point contact geome- 
try, which requires a much smaller gradient than 1/^s 
at the critical current density in the bulk, while the 
Josephson current is determined by the region within ^5 
from the junctions Since we choose the S and N lay- 
ers of the same transverse width, Ic/I^^^^ can be close 
to one for thin junctions. In such cases, current flow af- 
fects appreciably the superconducting order parameter 
[i.e., F{z) both inside and outside the N] and a self- 
consistent treatment becomes necessary (as is the gen- 
eral case pf| finding the critical current of a bulk super- 
conductoiES'Cj). Since for a clean SNS junction Rf^a? = 



[{2e'^/h){kj,/4Tr)]-^ « 1.58ha^/2e^ is just the Sharvin 
point contact resistance (i.e., inverse of the conductance, 
at half-filhng, per unit area a^ of our junction with infi- 
nite cross section), the absolute hmit of the characteris- 
tic voltage is /^""^-Rn — lAbAJe set by the bulk critical 
current /^""^ of the S leads ,t3 as shown in Fig. |l|. In 
three-dimensional (3D) junctions /^""^ = 1.09enA/hkF 
(per unit area a^, at half-filling) is slightly higher than 
the current density determined by the Landau depair- 
ing velocity Wdopair = A/hkp, at which supecfluid flow 
breaks the phase coherence of Cooper pairs,E3 because 
of the possibility of gapless superconductivity at super- 
fluid velocities slightly exceedingHj Wdcpair- Although our 
IcRn is always smaller than the ideal KO limit, it is still 
above the experimentally measured valueaJ in the inter- 
mediate junction thicknesses, which are about hundred 
times smaller than the KO limit. This suggests that addi- 
tional scattering confined to the interface regipB is indeed 
necessary to account for such small values.LJo 

Since self-consistent calculations require a phase gradi- 
ent inside the S (which we choose to be a boundary condi- 
tion in the bulk of the superconductor), we must carefully 
define how to parameterize the Josephson current. There 
are two possiMlities: either a global-phase change across 
the N regionc3 or the phase offsetEJ which is related to 
the phase change by a nontrivial scale transformation. 
We use a global phase change which in a discrete model 
like (|l|) requires a convention. The thickness of the junc- 
tion is defined to be the distance measured from the point 
zl, in the middle of the last S plane on the left (at zf ) 
and the first adjacent A'' plane (at z^ = zf + 1), to the 
middle point zn between the last N and first S plane 
on the right. Since 4'{z) is defined within the planes, we 
set (P{zl) = [0(^f ) + 0(-Zl^)]/2 to be the phase at zl, 
and equivalently for (jj^zn). The phase change across the 
barrier is then given by 



dz 



+ S(f){zii) - S(I){zl), 



(2) 



bulk 



where S(f>{z) is the "phase deviation" which develops self- 
consistently on top of the imposed linear background 
variation of the phase. The current versus phase change 
relation is plotted in Fig. |[ Non-self consistent jcalcu- 



lations predict that Ic occurs at 
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-K for bothQ ScS 



and long SNS junctions (at T — Q)u However, the self- 
consisteirt analysis leads to a sharp deviation from these 
notionsJi3 which is most conspicuous in our SNS geom- 
etry with a single normal plane. Moreover, even in the 
long junction limit (L = 60a ~ 15^o) we find 0c — 7r/2. 
The non-negligible A//i also leads to a lowering of (/)c (and 
a washing out of the discontinuities in /(0) at T = 0), but 
comparison with non-self consistent calculations, which 
take such normal scattering into account ,EJ shows that 
this is negligible compared to the impact of the self- 
consistency. 

The macroscopic phase of the order parameter 4>{z) 
varies monotonically (i.e., almost linearly) across the self- 
consistently modeled part of the junction. However, the 



plot of (5(/)(z), obtained after the linear background is sub- 
tracted from (/^(z), reveals a peculiar spatial distribution 
which depends on the thickness of the junction (Fig. ||). 
In the short and intermediate junction limits, 5(j){z) gives 
a negative contribution to 4'{z), which turns into a posi- 
tive one upon approaching I^. For thick enough junctions 
(e.g., L — 20a in Fig. ||) a small bump as the remnant 
of this behavior, persists at the SN boundary, but is 
completely washed out in the long junction limit. Thus, 
5(t){z) forms a "phase anti-dipole" (i.e., its spatial distri- 
bution has positive and negative parts opposite to that 
of the phase dipole, introduced in Ref. |l^), which is a 
self-consistent response to a supercurrent applied in the 
bulk. From the scaling feature of the phase-antidipole 
we conclude that such counterintuitive behavior of the 
"phase pile up" around the SN interface is generated by 
the finite A//i effects. 

Finally, we examine the local density of states (LDOS) 
p{ijj, Zi) on the central plane of the L = 10a junction, 
as shown in Fig. H. At zero Josephson current we find 
peaks in the LDOS, which are of finite width, corre- 
sponding to the Andreev bound statesc3 (ABS). More- 
over, instead of a non-zero LDOS all the way to the 
Fermi energy at w = (vanishing linearly as a; -^ 0), 
which stems from quasiparticles traveling almost paral- 
lel the SN boundary, a minigap Eg ^ A^//i is found 
which appears toJpe the consequence of finite A//i in- 
duced scattering.Ej The quantized bound states are the 
result of an electron (with energy below A) being retrore- 
flected into a hole at the SN interface, while a Cooper 
pair is injected into the superconductor, and vice versa. 
The hole is in turn transformed into an electron on the 
opposite surface, so that in the semiclassical picture, a 
bound state forms corresponding to a closed quasipar- 
ticle trajectory (i.e., an infinite loop of Andreev reflec- 
tions clectron^hole^electron . . . ) . The time-reversed 
ABS carries current in the opposite direction, and the 
two bound states are degenerate and decoupled (if there 
is no interface scattering). When the phase gradient is 
set within the S leads, a phase change appears across the 
junction (i.e., DC Josephson current), and the degenerate 
ABS split due to the Doppler shift. On the other hand, 
the minigap changes only slightly with increasing cj). The 
two Doppler split peaks drift apart monotonically until 
a bulk phase gradient corresponding to /c/2, when one 
of them reaches the BCS gap edge, while the other one 
approaches the minigap edge. The motion of the ABS for 
larger current be|G«mes more intricate and is described in 
detail elsewhere. Ea 
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